Quantum dynamics in light-harvesting complexes: Beyond the single-exciton limit 
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Primitive photosynthetic cells appear over three billion years prior to any other more complex 
life-forms, thus it is reasonable to assume that Nature has designed a photosynthetic mechanism 
using minimal resources but honed to perfection under the action of evolution. A number of different 
quantum models have been proposed to understand the high degree of efficient energy transport, 
most of them are limited to the scenario of single-exciton. Here we present a study on the dynamics 
in light-harvesting complexes beyond the single exciton limit, and show how this model describes 
the energy transfer in the Fenna-Matthew-Olson (FMO) complex. We find that the energy transfer 
efficiency above 90% under realistic conditions is achievable. 

PACS numbers: 05.60.Gg, 03.65.Yz, 03.67.-a 



I. INTRODUCTION 

Light-harvesting complexes consist of several chro- 
mophores mutually coupled by dipolar interactions re- 
siding within a protein scaffold. Due to their mutual 
coupling, light-induced excitations on individual chro- 
mophores (sites) can undergo transfer from site to site. 
Excitation energy transfer (EET) has been an interest- 
ing subject for decades, not only for its phenomenal ef- 
ficiency but also for its fundamental role in Nature [jj. 
Recent experiments on the exciton dynamics in photo- 
synthetic bio-molecules (for example, the purple bacteria 
and the Fenna-Matthew-Olson complex) have brought 
a long-standing question again into the scientific focus 
that whether nontrivial quantum coherence effects exist 
in natural biological systems under physiological condi- 
tions [2, 3]. In fact, evidence of quantum coherence has 
been found, suggesting that nontrivial quantum effects 
may be at the heart of its remarkable excitation trans- 
port efficiency 

Inspired by these experimental results, several stud- 
ies have attempted to unravel the precise role of quan- 
tum coherence in the EET of light-harvesting complexes 
[El-[l3j. and environmental decoherence and noise have 
been found to play a crucial role H3> EH ■ In these 
studies, the system is assumed to be initialized with a 
single excitation in site 1. This may not be realized 
precisely under experimental or natural operating con- 
ditions. Considering that primitive photosynthetic cells 
appeared over three billion years prior to any other more 
complex life-forms, it is not illogical to assume that na- 
ture has designed a photosynthetic mechanism using min- 
imal resources to gain maximal energy under the action 
of evolution. In this perspective, a model which allows 
a freedom to control the number of excitations in the 
complex at any time should be taken into account. 

This paper extends the theoretical formulation pre- 
sented in a recent paper [7[ to a scenario of multi-exciton 
and further examines issues relevant for realistic light- 
harvesting complexes. To this end, we identify the FMO 
complex with coupled cavities and introduce two fun- 
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FIG. 1: (A) The disordered structure of the Fenna-Matthews- 
Olson (FMO) complex. The FMO complex acts as an energy 
transfer channel in green sulphur bacteria guiding excitations 
from the lightharvesting antenna at site 1 to the reaction 
center at site 8. This picture is re-produced from [l5| . (B) 
The model to describe the FMO complex. Arrows between 
the cavities represent cavity-cavity couplings. Only couplings 
above 15 cm -1 are shown. The reaction center is modeled by 
the cavity, numbered 8, that is irreversibly coupled to the site 
3. It is worth bearing in mind that the coupling between cav- 
ity 1 and cavity 6 is week with respect to the other couplings 
shown in the figure. 



damental decoherence mechanisms (i.e., dephasing and 
dissipation) into the system. Under the semi-classical 
approach and with the quantum theory, we calculate re- 
spectively the excitation transfer efficiency. The deco- 
herence rates that are optimal to the ETE are found by 
numerical simulation of the equation of motion. With the 
optimal decoherence rates, the time evolution of popula- 
tion on each site is presented. We found that the optimal 
decoherence rates weakly depend on the initial excitation 
number on site 1. The non-local decoherence slightly al- 
ter the ETE, indicating that local decoherence dominates 
the energy transfer in the FMO complex. 
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The paper is organized as follows. In Sec. II we intro- 
duce the theoretical model for the FMO dynamics and 
define the excitation transfer efficiency (ETE) used here. 
Then, we analyze the ETE and the dynamics in the FMO 
complex with both semi-classical and quantum theories 
in Sec. III. Finally, we conclude our results in Sec. IV. 

II. MODEL DESCRIPTION 

In photosynthetic antennae, sunlight is absorbed by 
pigments and the excitation energy is transferred to the 
photosynthetic reaction center. The locations for these 
processes are physically and physiologically separated, 
suggesting that two-level systems are good to model the 
pigments with single-exciton but not enough for pigments 



with many excitations. To describe the light-harvesting 
process with many excitations on each site, we model 
the FMO complex by a coupled 8-cavity system, see Fig. 
[TJ The effective dynamics can be described by a Hamil- 
tonian describing the coherent exchange of excitations 
between chromophores or sites, 

7 7 

H = u i a \ a i + 9ij( a i a j + a * a ])i ! ^'> C 1 ) 

where (a,j) are the creation (annihilation) operators 
for site j, uij is the local site excitation energy, and c?y 
denotes the hopping rate of an excitation between the 
sites i and j. In the site basis, we follow [l6j and employ 
the Hamiltonian matrix elements (in units of cm -1 ) 
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Here the zero energy has been shifted by 12230 cm -1 for 
all sites, corresponding to a wavelength of ~ 800nm. Wc 
note that in units of % = 1, 1 ps _1 =5.3 cm -1 . Then by 
dividing <fe and uij by 5.3, all elements of the Hamilto- 
nian are rescaled in units of ps^ 1 . We can find from the 
Hamiltonian H that in the Fenna-Matthew-Olson com- 
plex (FMO), there are two dominating EET pathways: 
1 ->• 2 ->■ 3 and 6 -> (5, 7) -)• 4 ->• 3 (see figure [J). Al- 
though the nearest neighbor terms dominate the site to 
site coupling, significant hopping matrix elements exist 
between more distant sites. This indicates that coher- 
ent transport itself may not explain why the excitation 
energy transfer is so efficient. 

To obtain high energy transfer efficiency in the EET 
process, forward and backward energy transfer rates as 
well as the dissipation induced by the hot and wet sur- 
rounding environment need to satisfy a detailed balance 
condition. In the weak dissipation regime, the Lindblad 
master equation that is able to reliably describe exciton 
dissipative dynamics reads, 

-jj; = -i[H,p] + Cde P h (p) + Cdiss (p) 

+£-NLdeph(p) + C NL diss{p) + £ 8 (p) 1 ( 3 ) 

where Cdi S s(p) describes dissipation terms, and Cdeph(p) 
denotes dephasing terms. C NLdeph {p) and C NLdiss (p) 
represent non-local dephasing and dissipation, respec- 



tively. Here, 

7 

Cdiss{p) = Y [ ~ P a ] a 3 ~ a ) a jP) > ( 4 ) 

7 

Cdeph(p) = } j jj(2rijpnj - prijUj - njnjp), (5) 
i=i 

7 

£-NLdiss(p)= Y r ij(2a l/ 9a]- / oa]a i -a]a i p),i ^ j, (6) 
and 

7 

jO-NLdeph(p) = Y lij( 2n iP n j ~ P n J n i ~ njflip),i ^ j, 

(7) 

where rij = a]flj denotes the exciton number operator 
for site j. 

The local decoherence may come from the couplings 
of the sites to individual environments, whereas the non- 
local decoherence CNLdeph(p) and CNLdiss{p) can be un- 
derstood as a result of the interaction between the site 
and a common environment. From a quantum physical 
perspective, environmental radiations whose wavelength 
is larger than the length-scale of the FMO molecule (~ 8 
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FIG. 2: Energy transfer efficiency versus F% without decoher- 
ence (Fj = fj = Tij = -y i3 = 0, i, j=l,2,...,7). (a) For semi- 
classical approach and, (b) for quantum theory. iVo = 100 is 
taken for this plot. 



nm) are reasonably assumed to take the role of a com- 
mon environment. Thus the decoherence considered here 
is reasonable. 

Recent work [3| suggests that it is the site 3 that cou- 
ples to the reaction center. The total transfer of excita- 
tions into the reaction center is measured by the popu- 
lation in the center, numbered 8, which is populated by 
an irreversible decay process with rate r 8 from the site 
3. We phenomenologically model this irreversible process 
by the Lindblad operator, 



£s(p) = r 8 (24a 3j o4c 



pa\a 3 a s a\ 



ala 3 a 8 alp) . (8) 



To match the observation in experiments, we assume 
there are No excitations initially in site 1. The model is 
completed by introducing a quantity by which we mea- 
sure the energy transfer efficiency. The rescaled popula- 
tion in the reaction center given by 



n 8 (T) _ Tr(a|a 8 p(T)) 



Nq 



N 



(9) 



at a specific time T is good for this purpose. 

Note that the reaction center is directly connected only 
to site 3. The decoherence rate r 8 plays an essential role 
in the excitation transfer. In the next section, we will use 
the decoherence rates to optimize the transfer efficiency 
defined in Eq.©. Several cases are considered, in each 
case Tg can not be zero, because null T 8 leads to zero 
transfer efficiency. 



ered. In the semi-classical approach, we will use the ap- 
proximation (<zjaj) = \a-j\ 2 with ctj = Tr(paj). Whereas 
in the quantum regime, we approximate (a\asala 3 ) as 
(a 8 a s) (o^ 3 a 3 ) . By these approximations, we can derive 
the equation of motion for the system and calculate the 
energy transfer efficiency. We first focus on the case 
where only local decoherence exists, namely Fjj = and 
7y = 0, then move to the case with non-local decoher- 
ence. The energy transfer time is taken to T = 5 ps _1 , 
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FIG. 3: (color online) Efficiency as a function of time, only 
local decoherence is considered. Figure (a) is the results with 
semi-classical approach, where black solid line is for the null 
decoherence case (Fj = 7j =0), and blue circle and red 
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respectively, j = 1,2,..., 7. Other parameters chosen are 
Fg — 1.94 and iVo = 100. (b) is the efficiency from quan- 
tum theory. Blue solid line stands for optimal efficiency with 
rates (71,72,73,74,75,76,77) =(0.74, 24, 0, 5.2, 50.6, 0, 15), 
Fg = 0.32 and the other Fj = 0.0005. The efficiency reaches 
its maximum 91.77% that is optimized at 5 ps. Red circle 
and black triangle lines represent the null decoherence case 
with Fs = 0.9 (the corresponding efficiency is 67.06%) and 
Fg = 0.32 (corresponding to efficiency 60.73%), respectively. 



which is the relevant time scale in recent experiment [__J . 
Our results suggest that it is the careful interplay of quan- 
tum mechanical features and the unavoidable couplings 
to environment that will leads to the optimal system per- 
formance. In particular, the local decoherence dominates 
and actually helps the excitation transfer, as we will show 
below. 



III. ENERGY TRANSFER EFFICIENCY 



A. Semi-classical approach 



Having this model, we now explore the energy trans- 
port in the FMO complex. Two approaches are consid- 



Define ctj — Tr(paj) and rig = Tr(aga 8 p), the master 
equation Eq. §5§ yields, 
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n 8 = 2r 8 |a 3 | 2 (n 8 + 1), 



where the dot represents time derivative, and the time ar- 
guments have been omitted to shorten the representation. 
To obtain Eq. (|10[) . the approximation (aja.,) = \otj\ 2 in 
the last equation of Eq. (flOl) has been made. These 
equations compose a closed set of equations governing 
the dynamics of the FMO complexes. 

As shown in [7], a completely coherent dynamics is 
often not most ideal for the excitation transfer from 
the chromophores to the reaction center, and hence the 
coherence solely can not explain the very high exciton 
transfer efficiency observed in experiments. Our model 
confirms this result. In fact, with Tj — jj — Fy — jij — 
= 1,2,..., 7), the transfer efficiency arrives at its 
maximum 0.625 with Tg = 1.94 ps _1 (see Figs. HJa) and 
Ha)). 

Further numerical simulations show that under the 
semi-classical approximation, dissipation and dephasing 
play similar roles in the dynamics of FMO, see Eq. (|TU| . 
Numerical simulations show that neither dissipation nor 
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where a = 1 for m — 3 or n — 3, and a — 2 for m = n = 3, 
otherwise a = 0. Eq. (fTTjl is accurate when (ogOgaJaa) 
= (a s as) (alas) . In fact, our numerical simulations show 
that this is exactly the case for small number of exciton. 
For instance, with N — 1 or 2, Monte Carlo simula- 
tions show that (ogOgOgOs) — (alas) (0303) = 0. On the 
other hand, our model backs to the two-level model for 
the FMO complex when N = 10, and the numerical 
results given by Eq. (jlip is in agreement with that in 
0. It is difficult to prove (0,^0303) — (agOs) (03(13) = 
for any Nq. Fortunately, the site 3 and the reaction cen- 
ter 8 are connected only through the irreversible process 
C&(p), and this term unlikely creates entanglement be- 
tween site 3 and 8. Therefore if site 3 and 8 are ini- 
tially in a separable state, then they will remain unen- 
tangled forever. For site 3 and 8 in separable states, 
(alasalas) — (ogOg) (0303) holds true. This is not a 
proof, so we prefer to treat (ogOgOgOs) = (ctgag) (03(23) 
as an approximation. 



(10) 

I 

dephasing can improve the energy transfer efficiency un- 
der the semi-classical approximation (see Fig. [3](a)). 
This is due to the fact that dephasing and dissipation 
play a same role in the dynamics under the semi-classical 
approach, see Eq. (fTU|) . Remind that the local dissipa- 
tion always spoil the exciton transfer from sites to the 
reaction center, it is not difficult to understand why the 
decoherence can not help the exciton transfer within the 
semi-classical approach. 



B. Quantum theory 

We now examine whether the quantum theory can ex- 
plain the high exciton transfer efficiency. For this pur- 
pose, we approximate (ogOgaJos) as (agdg) (030,3). Define 
n mm = TrQtd^am), and n mn = Tr(p<4ja„). We obtain 
fromEq. © 



(11) 



We first study the case where only local decoherence 
exists, namely, = 7.^ = 0. By numerically solving 
Eq. (jllj) . we find that without decoherence the maxi- 
mal energy transfer efficiency is 67.06% (see Fig. [2fb)), 
which is a bit larger than that given by the semi-classical 
approach. However, with the assistance of local decoher- 
ence, an efficiency over ninety percent can be obtained 
(See Fig. EKb)). In order to obtain such a high energy 
transfer efficiency, we have to optimize the 15 decoher- 
ence rates (Tj, 7^, and Tg, j — 1, 2, 3, 7), this is a 
time-consuming task. To save the computer time, we 
specify the dissipation rates to Tj = 0.0005 for each site 
to match the measured life time in experiments, which is 
the order of nanosecond [3j for the light-harvesting sys- 
tem. The optimal dephasing rates for highest efficiency 
are found numerically as (71, 72, 73, 74, 75, 76, 77) =(0.74, 
24, 0, 5.2, 50.6, 0, 15) and T 8 = 0.32, the corresponding 
efficiency is 91.77% for initial Nq = 100 on set l.(see Fig. 
|3]). Further numerical simulations show that the energy 
transfer efficiency can reach over ninety percent for al- 
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FIG. 4: Population ratio of j-th site defined by Pj = as a 
function of time (number 1,2,. ..,7 label the sites), (a) is for the 
semi-classical approach. And (b) is for the quantum theory. 
Dissipation and dephasing rates are the same as those which 
optimize the efficiency with initial exciton number No = 100. 




FIG. 5: (color online) Efficiency versus excitation number TVo 
Parameters are chosen in such a way that the efficiency is op- 
timized with No — 100. (71,72,73,74,75,76,77) =(0.74, 24, 
0, 5.2, 50.6, 0, 15), F 8 = 0.32 and T = 0.0005 is for the red 
triangle line. In contrast, we plot the results without deco- 
herence in green square and blue circle lines with F$ = 0.32 
(quantum) and T$ = 1.94 (semi-classical), respectively. The 
later two Tg optimize the transfer efficiency for null decoher- 
ence case. 
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most arbitrary exciton number Nq (from 1 to 10000) on 
site 1, though the optimal dephasing rates are different. 
With the optimal dephasing rates, the ratio of exciton 
number on each site to the total number of exciton (de- 
fined as p m — <s ^f L ) is plotted, see FigfJ] (b). We find 
that after a rapid increasing (except for site 1), the exci- 
ton number decreases for each site, indicating that most 
of the exciton are transferred to the reaction center at 
site 8 (see Fig. S{b)). In contrast, the dynamics from 
the semi-classical approach (see Fig. H(a)) is also plot- 
ted. More oscillations can be observed, leading to less 
population transfer to the reaction center. 

It is interesting to note that the optimal decoherence 
rates obtained for different exciton number Nq are ap- 
proximately the same values. To show this point we have 
plotted in Fig. [5] the efficiency as a function of the exci- 
ton number iVo, the decoherence rates used in this figure 
are optimal for No = 100. Clearly, the decoherence rates 
optimal for No = 100 can also result in high transfer 
efficiency for a wide range of No. For instance, the effi- 
ciency remains above ninety percent at No = 200 with 
the optimal decoherence rates for Nq = 100. This obser- 
vation suggests that the high efficient energy transport, 
which has been found for a fixed excitation number re- 
mains valid for other experimental and natural operating 
conditions. Biologically, this means the light-harvesting 
system is robust against the number of photon captured 



FIG. 6: (color online) Efficiency as a function of local dis- 
sipation rate =T,j — 1,2,..., 7). This figure is plotted 
to show the effect of non-local decoherence on the transfer 
efficiency. Red solid line represents a case with nonlocal dis- 
sipation (Fij = r), whereas the blue dash line is for the case 
without nonlocal dissipation (Ty = 0), i,j = 1, 2, 7. Other 
parameters chosen are Ps = 6, jj = 0, and 7^ — 0. 



at site 1. Meanwhile, we observe that the dependence on 
small No (see Fig. [5J is stronger than the dependence on 
large Nq, suggesting that the transfer efficiency is sensi- 
tive to the variations of small No. This feature can be 
understood by carefully examining Eq. (jlljl . Clearly, the 
losing rate of population on site 3 depends on the popu- 
lation in the reaction center, n^- As a back action, the 
gaining rate of population on the reaction center depends 
on n 33 . 



C. Effect of non-local decoherence 

In this section, we examine the effect of non-local deco- 
herence on the energy transfer efficiency and the dynam- 
ics of the FMO complex. When taking Eqs. ([BJ and (O 
into account, we must guarantee that the non-local deco- 
herence terms with rates and 7^ should keep the pos- 
itivity and trace preserving of the density matrix Q ■ For 
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this reason, we choose all Tij and 7^ positive to optimize 
the transfer efficiency. We find that the energy transfer 
efficiency can be increased by properly choosing non-local 
dephasing rates. For example, with the optimized local 
decoherence rates ({71,72,73,74,75,76,77} ={0.74, 24, 
0, 5.2, 50.6, 0, 15}, T 8 = 0.32, and T = 0.0005), the ef- 
ficiency can be increased from 91.77% without non-local 
decoherence to 91.907% with non-local decoherence rates, 
(7i7,77i,725,75 2 )=(0.74,0.74,24,24) and 7y = 0,{i,j} + 
{1,2,5,7}. Extensive numerical simulations show that 
energy transfer efficiency can not be increased by weak 
non-local decoherence (orders of ns _1 ), but strong 
non-local decoherence can improve the energy transfer 
efficiency as Fig. |6] shows. 

IV. CONCLUSION 

In summary, we have studied the dynamics in light- 
harvesting complexes beyond the single exciton limit and 
optimized the energy transfer efficiency in the Fenna- 
Matthew-Olson (FMO) complex. To describe the multi- 
excitation scenario, we have proposed a new model for the 
propagation of excitation transfer, the model consists of 



7 coupled cavities and a reaction center. Four types of 
decoherence, including local dephasing, local dissipation, 
nonlocal dephasing and nonlocal dissipation are consid- 
ered. To match the life-time of exciton observed in exper- 
iment, we fixed the local dissipation rates to 0.0005 ps _1 . 
The local dephasing rates that optimizes the transfer ef- 
ficiency are given by numerical simulation. We find that 
for multi-excitation case, the energy transfer efficiency 
can be over 90% under realistic conditions. Non-local de- 
coherence can slightly increase the efficiency, but it seems 
not important in the light-harvesting mechanism. More- 
over, we find that the transfer efficiency is not sensitive 
to the initial excitation number at site 1. This suggests 
that the light-harvesting antenna may capture more pho- 
tons once and the experimental conditions are flexible to 
simulate the light-harvesting in FMO complex. 
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